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Abstract 

We examine the neutrino mass matrix in the version of Zee model where both 
Higgs doublets couple to the leptons. We show that in this case one can accommo- 
date the large mixing angle (LMA) MSW solution of the solar neutrino problem, 
while avoiding maximal solar mixing and conflicts with constraints on lepton fam- 
ily number- violating interactions. In the simplified scenario we consider, we have 



the neutrino mass spectrum characterized by mi ~ m2 — y Am^^j^/ sin 29 and 
m^/m\ ~ cos 26, where 6 is the solar mixing angle. 
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1 Introduction 



The results of the atmospheric PP and solar neutrino experiments (for a recent presentation 
of the solar results see, e.g., Ref. |2]) indicate that neutrinos have non-zero masses. A 
suitable environment for obtaining Majorana neutrino masses is to extend the Higgs sector 
of the Standard Model ^3J. Extending the Standard Model with a singly charged gauge 
singlet scalar and adding a second Higgs doublet allows to write down an explicitly lepton 
number-violating interaction in the Higgs potential and leads to 1-loop neutrino masses 
0]. In the following we will only consider the lepton sector of the Zee model. 

The Zee model is traded in two versions in the literature: the original model |3], 
which we will call general Zee model (GZM) in this paper, and a simpler version of 
the Zee model where only one of the two Higgs doublets couples in the lepton sector 
1^. The latter version, which we will call restricted Zee model (RZM), can naturally 
be achieved with a discrete symmetry |6j it has the advantage that the family lepton 
number-violating interactions mediated by the couplings of the three physical neutral 
scalars are absent. The interesting point is that the RZM leads to a symmetric neutrino 
mass matrix with zeros in the diagonal^; this mass matrix is called Zee mass matrix in the 
literature. One can easily check that all phases in the Zee mass matrix can be absorbed 
into the left-handed neutrino fields and, therefore, in neutrino oscillations no CP violation 
is observable if the Zee mass matrix is the correct neutrino mass matrix. 

The Zee mass matrix has a special feature: it has been shown [3 |H] that it allows to 
accommodate only bimaximal mixing, i.e., solar and atmospheric mixing angles are both 
very close to 45°. While this is in excellent agreement with the atmospheric neutrino 
data, results from the solar neutrino experiments are not that well compatible with a 
solar mixing angle 9 ~ 45°, though this value is also not excluded [21101 • The purpose of 
the present paper is to show that in the GZM it is no problem to accommodate the large 
mixing angle (LMA) MSW solution of the solar neutrino deficit, while at the same time 
all constraints on the additional couplings in the lepton sector stemming from the second 
Higgs doublet are respected. In the GZM also the diagonal elements of the neutrino mass 
matrix are non-zero in general. Since the GZM is a quite rich and intricate model, we 
restrict ourselves rather to an "existence proof" of the LMA MSW solution within the 
GZM instead of discussing the GZM in full generality. 

Let us outline our procedure: 

Point 1: We assume that Ues = (t/ is the neutrino mixing matrix) and the atmospheric 
mixing angle is exactly 45° ^U] . Then in the mixing matrix the only free parameter 
is the solar mixing angle 6. This gives a certain form of the mass matrix with 4 
complex parameters (Section EJ. 

Point 2: The general problem is still quite intricate, so we set one of these parameters of 
the mass matrix equal to zero and assume that the remaining ones are real. In this 
scenario we can relate in a simple way the three real parameters with the physical 
quantities 6, the solar mixing angle, AtjIq = \m1 — m^l, the solar mass-squared 
difference, Am1^^ = jmg — (m^ +m2)/2|, the atmospheric mass-squared difference 

^We always work in a basis where the charged lepton mass matrix is diagonal. 
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and m3 {rrij with j = 1,2,3 are the neutrino masses) and it is possible to have the 
LMA angle solution of the solar neutrino problem^ (Section E)). 

Point 3: Now the GZM is brought into play. After a discussion of the neutrino mass 
matrix in this model (Section 0} we assume that all quantities in the model are 
real and we set to zero all but two of the additional Yukawa couplings present in 
the GZM; this has the purpose of avoiding as many family lepton number-violating 
neutral scalar interactions a possible. We show that with these two additional 
coupling constants the restricted mass matrix mentioned in Point El which allows 
for the LMA MSW solution, can be accommodated (Section Hj). 

Point 4-' We complete our procedure by a numerical discussion of the parameters of our 
scenario and by estimates of the rates of (family) lepton number-violating processes 
(Section ini). 

We also review the features of the RZM as a limit of the GZM (Section Ej) and summarize 
the results (Section |H)). In the appendix we present the general formulas for the 1-loop 
Majorana neutrino mass matrix induced by charged scalar loops. 



2 Neutrino mixing and the mass matrix 

The Majorana neutrino mass matrix Ai^ is diagonalized with a unitary matrix V by 



V A^j^V = m = diag (mi, m2, ma) . (1) 

With the assumptions mentioned in the introduction in Point ^ we can write the matrix 
V as 



V = e'^Ue''^ with U 



/ cos 9 sin 6* 

— sin 9/y/2 cose/V2 
V sme/V2 -cose/V2 1/V2 



The phases in V suggest the definitions 



M'^ = UflU^ with Ml = e'^M^e'^ and = e 



-2i0 



m . 



(2) 



(3) 



Then Ai'^, can be expressed by fi and the parameters of U in the following way: 
M' 



-Csifli - IJ,2)/V2 (sVl + cV2 + /i3)/2 (-sVl - + /i3)/2 

V Cs(/ii - /i2)/V2 (-sVl - + /i3)/2 (sVl + cV2 + /i3)/2 



(4) 



with c = cos^ and s = sinO (see, e.g., Refs. [HE])- Thus, the mass matrix has the 
structure 

a —b b 

M:,= \ -b c d \ . (5) 
b d c 



^This means that 9 is large but safely smaller than 45°. 
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Consequently, for having 11^^ = and an atmospheric mixing angle of exactly 45° the 
following conditions on Ai'^^ are necessary: 



Condition 1: M'^^^ + M'^ 



0, 



Condition 2: M', 



V llfj, 



M' 



Conditions: tan2e = 2v^- 



e R. 



a — c + d 

Finally, with the parameterization the complex masses fij are found as 



a + c-d± {a-c + df + 8b'^ 



2 

c + d . 



a + c - dT {a - c + dY + 8b'^ 



1/2 



1/2 



(6a) 
(6b) 

(6c) 



(7a) 

(7b) 
(7c) 



3 A simplified mass matrix 

Having discussed the general form of the mass matrix which leads to the mixing matrix 
((21), we now investigate the consequences of the following simplifying assumptions: 

a,b,d eH and c = . (8) 



In the next section, this scenario will be reproduced in the framework of the Zee model. 
With the reality assumptions, the quantities fij ^ are identical with the neutrinos masses 
apart from possible signs. The experimentally accessible quantities are expressed as 

b 

tan 20 = 2^2 -, (9a) 

a + d 

r n 1/2 

Am| = \a-d\ {a + df + Sb'^ , (9b) 

Am^,^ = ]^{a' - d^) + 2b\ (9c) 

ma = \d\ (9d) 

by the parameters a, 6, d. We have chosen as representative of the absolute neutrino 
mass values, since it is simply given by Eq. ((9H|l . Without loss of generality we will adopt 
henceforth the following conventions: 0° < < 90°, mi < m2, /is = m^. It follows from 
the last relation and from Eq. (|7c)) that d is positive. Note that in Eq. no absolute 
value of the right-hand side of the equation is necessary, because it must be positive. The 
argument goes as follows. Suppose that (a^ — (i^)/2 + 26^ = — Am^^j^. Then it follows that 
d^ — > 2Am^tjjj. Therefore, d^ — a? is positive, which allows to derive the inequality 
Am| >d^-a^> 2Am^^^ from Eq. (jHEI). This is a contradiction to the values of the 
mass-squared differences, fitted from the data |2j. 

In Eqs. (jHij), (j9H)), and (j9H|), four physical quantities are expressed by three 
parameters. Therefore, a consistency condition exists, which is given by 

Am^tj^ = ^r/Amll cos 2611 + ^ {^\Jml + ri^m%\ cos2d\ + rj'm^ tan^ 26* . (10) 
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The signs t] and r]' occurring in this equation are rj = sign (a^ — c?^) and rj' = sign a. In 
the context of the Zee model we will finally need the relations 



= m3 + r^AmQl cos26'| , (11a) 
= ^Aml^-^r]Aml\cos29\. (lib) 

Looking at the consistency condition (|lUj) and assuming that is of the order of 
Amgl cos 26*1 or smaller, we obtain sin^ 2^/| cos20| ~ Aml^^/ArriQ, which amounts to 
bimaximal mixing. Since we want to show that the Zee model allows to avoid bimaximal 
mixing we concentrate on 

> Am|| cos26'| . (12) 

With this assumption it is easy to obtain an approximate expression for mg. One can 
check that for r]' = —1 one arrives again at bimaximal mixing. Using r/' = 1, we have 
a > and we easily calculate 

2.2(1 1 Ami |cos20| 1 / Ami V o ^ , 

ml = Am^_ <^ — - 7]-- — — 5 — ^ + — - — ^ s n^ 29 + ...)■ . 13 

^ ^*"'|tan2 2^ '2 Am^,^ sin2 2^ 16 V Am^j^y ' ^ ^ 



Using this equation the condition (|T^ implies 



cos2^|»,— ^. (14) 



.0 

Aml^^^ 



Equation ()13|) together with a > and ci > 0, allows to estimate from Eq. ()lla|) that 

, 1 Ami, ^ , , 

a — a ^ -7] ^ = sin 26^. (15) 



2 . //^^2 



atm 

From Eqs. ()9a|) and ()9b|) and the convention mi < m2, we can express the masses mi and 
m2 as 

1 / Am^ \ 

Inspection of Eqs. ()7a|l and ()7bjl reveals that our conventions fix the signs of fii^2- sign /ii = 
— sign/i2 = —T]. Then, with Eq. p5|l . an estimate of the neutrino masses which neglects 
the solar mass-squared difference is given by 



Am^t^ ^jAml^ 
mi ~ m2 ^ . and ma ~ -j- — — — . (17) 

sm26' I tan 26' I 

This equation tells us that m^, < mi 2, at least in the regime of large solar mixing. 
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4 Neutrino masses in the general Zee model 



In the previous section we have discussed the mass matrix determined by Eqs. (0) and 
(jS)) without reference to any specific model of neutrino masses. Now we introduce the Zee 
model |S] and discuss the neutrino mass matrix in the case that both scalar doublets of 
the Zee model couple in the lepton section. The Yukawa Lagrangian is given by 



Cy = -Y. LVact^dR + L^C-HT2fLh+ + h.c, 



a=l 



where / is an antisymmetric 3x3 matrix 1,4,. The mass matrix of the charged leptons 
arises at tree level through 





1 



and Mi> 



(19) 



with V = Y|fiP + |f2P — 246 GeV. The physical charged scalar fields H^, H2 with 
masses Mi, M2, respectively, and the would-be Goldstone boson 0+ are obtained by the 
unitary transformation IHl ^| 



h- 



112 !hy^ !hiY_ 



12 



V 







w, 



21 



w, 



22 




(20) 



As anticipated in the introduction, we assume to be in a basis where the charged lepton 
mass matrix is diagonal, i.e., = M^, where the hat symbolizes that this mass matrix 
is diagonal. In this basis we have 



Ti = —{V2Mi-V2T2 



(21) 



In our parameterization the fiavour-changing Higgs couplings are given by the off-diagonal 
elements of r2. Furthermore, off-diagonal elements of the 2x2 unitary matrix W are 
present because of the vacuum expectation value of the lepton number-violating term 



Xh 



'2 + h.c. 



(22) 



in the Higgs potential. 

The physical charged Higgses couple to the leptons in the following way: 



ViV 



Vl 



H^Ln2f)lL{W2lHt + W22Ht) + h.c. 

With the formulas in the appendix we obtain 0] 

2 . . . -t 

Mu = Y. A^MJiMf, Mf)A^j^ + transp. 



(23) 



(24) 
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with 



= (~V2—Me + —T2) Wii and A{ = 2fW2i ■ (25) 

\ ViV Vi J ' 

The infinity in Eq. cancels |lj because ELi W2jW^j = 0. Defining 



from Eq. ()24|1 we obtain the final result jH EI 

X I v^-^ (Ml J/ - /JM|) - ^ (7M, J/ - / JM,7^) I (27) 

with J = diag (J(Mf,M|,m2)) (a = e,/x,r) and 7 = T^. Note that for < Ml^ the 
function J simplifies to 

J~ 21n(Mi/M2). (28) 

For the product W21W11 of elements of the charged-scalar mixing matrix W, we obtain 
the relation 

W21W* = ^ V T- ■ (29) 

y2(Mf-M|) ^ ^ 

It shows explicitly that the Majorana neutrino masses are proportional to the coupling A 
in the Higgs potential (j22I)- 

In the GZM considered here, there are family lepton number-violating processes in- 
duced by the charged and the neutral scalar interactions. Experimental bounds constrain 
the coupling matrices / and 



5 The simplified mass matrix within the Zee model 

In this section, our aim is to reproduce the neutrino mass matrix defined by Eqs. (0) and 
(jHl) within the GZM. In order to save the amount of writing we introduce the notation 

M-iy = A (r^ f — ff^^ — B (^"yff — fr'y'^^ with f = diag (me, m^, m^)/T; . (30) 

Both constants A and B are of the order of 1 GeV, resulting from dividing the electroweak 
scale by IGvr^: 

B = 2W2iW;,x——\n-^x- and A = V2^B. (31) 
(47r)^ M2 vl V 

Since the off-diagonal elements of 7 introduce fiavour-changing neutral interactions, we 
adopt the philosophy to set to zero as many of them as possible. As we will see, it turns 
out that 

7e. 7^0 and 7rr7^0, (32) 
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and all other elements of 7 being equal to zero, is sufficiently general to avoid bimaximal 
mixing, which necessarily happens for 7 = ||7j. It can easily be checked that with this 
assumption Condition 2 (jFB]) is fulffiled by having Ai^^^ = M-'utt — 0- Furthermore, we 
assume that all quantities we deal with are real: A, B, the elements of the matrix /, and 
7er and ■jrr- Thus we identify M.^, with M.^. 

We have to reproduce with the mass matrix the parameters a ()lla|l . b (jllbj) . d 
(|9H]) and Condition 1 (jUaf) . Thus we have five coupling constants, three in / and 7er, 
■jrr and four relations. It is convenient to express the other four coupling constants as a 
function of •jrT- 

f - fi I A{r^^-rl) - Brr-frr \ . . 

^(^2 _ ^2) I ^ + A(r2 - r2) - i?r,7.. J ' ^ ' 

= A{tI - rl) - BVrlrr ' ^^^^^ 



ITT 

a 



7. 



er 



2Brrh 



{A{rl-rl)-Brr^,,] . (33d) 



6 Numerical estimates 

Let us now estimate the values of the coupling constants. For definiteness we take 
^"^atm = 3x 10~^ eV^ and A = 2 GeV. Furthermore, we need the values of ~ 5.22 x 10"^ 
and mfj^/m'^ ~ 3.535 x 10"'^. Defining x = B'^-rr/Arr and assuming that^ x ~ 1 and 
(1 — mil ml — x)/{l — rn?^/m1 — x) ~ 1, we obtain 

^ -3.7xl0-^x^, (34b) 
\ — X 

f,r ^ -5.2x 10-^x } -, (34c) 

I tan 26\ (1 — xj 

7,. ^ -7.2x lO-^x -. 34d 

V tan2t^ 

These equations serve to see the orders of magnitude and any effects of AttIq are neglected. 
As can be seen from Eq. fll5|) . a considerable amount of finetuning is involved in order to 
reproduce the solar mass-squared difference. 

Now we concentrate on the LMA MSW solution of the solar neutrino problem, where 
6 is in the ffist octant. With Eq. ()9a|l it follows that 6 > 0. In this case a representative 
value of the mixing angle 9 is given by the best ffi value sin^2^ ~ 0.75 (tan 26' ~ 1.71) 
of Ref. with the corresponding mass-squared difference AttIq ~ 3.2 x 10~^ eV^. We 
note that in this case we have |/er//e/i| ~ ''^^/'^r; which is similar to the case of the Zee 
mass matrix [7j. On the other hand, in our scenario we have /er ~ //^r; whereas for the 

■^This avoids some finetuning for Jtt- Note, however, that a; = is also possible. 
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Zee mass matrix the relation |/e^| ^ |/^^| holds [7j. As far as 7 is concerned, with x ~ 1 
we can have ~ 7rT 10~^ ^ 10^^. 

Due to our assumption (jH^ . flavour- changing neutral scalar interactions at the tree 
level are very constrained. Among the charged lepton decays we only have r~ — » e^£^£^ 
with i = e, fi. A generous estimate of the branching ratio of this decay for i = fi is 
obtained by 

B.(.-^e-.V)^g^.i(^)%10-, (35) 

where we have taken G^^Mq ~ 10~^ (Mq ~ 100 GeV) and Mq is a generic neutral Higgs 
mass. At the 1-loop level, neutral Higgs exchange also induces the decay e~7. 
Making again an estimate, we obtain ^Hl 

B.(.-^e-,)^J-^(-,„,„)=<10-. (36) 

Both estimates are well compatible with the experimental upper bounds on these branch- 
ing ratios of the order of 10~^ llSj. More detailed discussions of these decays are found 
in Refs. OUn]. We have used |7e^| < lO^^ and |7^^| < lO'^ in Eqs. §^ and 

Also the charged scalars participate in various charged lepton decays as intermediate 
particles. Numerous decays of the type ^ £2 + 2 neutrinos proceed at tree level via 
the Lagrangian ()23|1 . According to the couplings in this Lagrangian we can distinguish 
between /, 7 and f vertices, and we can have decay amplitudes with all possible combi- 
nations of these vertices, except 7-7, which is forbidden due to the restricted form of 7, 
Eq. (jH^ . E.g., to the Standard Model amplitude of ordinary muon decay there is an am- 
plitude with couplings r^rg and another one with ferffir] another example is r~ u^i^^fi' 
with couplings ferf^^i, which is lepton number violating. The branching ratios of all these 
decays are negligible because of the smallness of the coupling constants fa/3 and ■jer and 
the ratios Tq,. Also neghgible are radiative decays — * ^2 +7 induced by charged Higgs 
loops jni- In this case one always has two /-vertices in the loop graph, except in the 
case of the amplitude for ^ e~ + 7, where there are two contributions, proportional 
to fenf^iT and 7er^e- Reccut reviews of the restrictions on the coupling constants fa/s are 
found in Refs. 

Scalar contributions to the anomalous magnetic moments of the electron and muon 
involving / |THj and f-couplings are totally negligible because these constants are too 
small. A contribution from the 7-couplings to the electron magnetic moment coming 
from T exchange is proportional to 'yerlre il9j and is thus zero in view of 'fre = 0. 

In our scenario we have Aiuee = d 7^ 0. The matrix element Aiuee is identical with 
the effective neutrino mass {m^) probed in neutrinoless double-beta decay. Therefore, 
this decay is allowed and the effective neutrino mass is given by 

0.05 eV , , 

This represents an order of magnitude which is accessible in future experiments (for recent 
reviews see, e.g., Ref. [20] )• 
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Having seen that the LMA MSW solution of the solar neutrino problem can be accom- 
modated in our scenario, we now proceed to the small mixing angle (SMA) MSW solu- 
tion. In this case we take as illustration the best fit value of Ref. [2], sin^2^ ~ 2.3 x 10~^ 
(tan26' ~ 4.8 x 10^^), which has a corresponding Am| ~ 0.5 x 10^^ eV^. From Eq. ()34a|) 
we see that now /e^ ~ 0.05 becomes relatively large and barely compatible with the re- 
quirement that the Zee boson does not have an effect on the muon decay rate so that it 
does not destroy the agreement in electroweak precision tests jTTj. Thus in our simple 
scenario we cannot incorporate safely the SMA solution. 



7 The limit 







To make contact with Refs. [H E], we explore the effect of r2 = 0. This means that 
all diagonal elements of the neutrino mass matrix are zero j3] and, in our notation, we 



have a = 0, i.e., 



cos2^^| and rj = —1 (see (Eq. (jllap ). These relations and 



inspection of the consistency condition (jlOp leads to 



7713 



1 

~ 2 



Am| 



rrii ~ m2 



and I cos 26* I 



1 Am| 
4 Am2,^ 



(3^ 



From the last relation we read off that 9 is 45° for all practical purposes. The consequences 
for the coupling matrix / are 



fe 



m 



m 



and 



fe 



f> 



(39) 



The first of these two relations is obtained by taking the ratio of Eq. (j33bjl and Eq. ()33aj) . 
where a = is taken into account. In the second relation, Eq. ()llb|) has been used. The 
results dSHl) and ^ agree with those of Refs. [HIHl. 



8 Summary 

In this paper we have discussed neutrino masses and mixing in the general Zee model, 
where both Higgs doublets couple in the lepton sector. In this endeavour, we were mo- 
tivated by the result that for r2 = 0, where r2 denotes the Yukawa coupling matrix of 
the second Higgs doublet, the Zee mass matrix leads to bimaximal mixing. It is true 
that a mixing angle of 45° is perfect for the description of the atmospheric neutrino data, 
but it does not represent a very good fit for the solar neutrino data. The general Zee 
model is a rather rich and intricate model. Therefore, in order to simplify the analysis, 
we have assumed that all quantities appearing in the neutrino mass matrix are real and 
that in the second Yukawa coupling matrix r2, which is set to zero usually, only two 
elements, the er and rr elements, are non-zero. We have shown that this is sufficient to 
accommodate solutions of the solar neutrino problem with a large mixing angle instead of 
maximal mixing; on the other hand, in the atmospheric sector maximal mixing remains. 
At the same time, all dangerous processes induced by the new couplings are sufficiently 
suppressed. Actually, we could have even set the rr element of r2 equal to zero and used 
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a single non-zero element in this coupling matrix. However, the rr element represents a 
possibility to achieve a further suppression of the potentially more dangerous er element. 
As in the case r2 = 0, the neutrino mass is the smallest neutrino mass, and the result- 
ing mass spectrum is rather of the type which is called "inverted hierarchy". However, 
whereas for r2 = one has <^ nii ~ m2, now ~ y^Am^tm/ tan 2^ is of the same 
order of magnitude as the other two masses. 

Numerous finetunings are involved in our scenario. Implementing the condition (j6a|) 
and setting all but two elements of r2 equal to zero represent rather severe finetunings. 
These procedures were useful for exactly having Ues = and an atmospheric mixing angle 
of 45°, and for avoiding a dangerous class of lepton family number- violating interactions. 
In order to reproduce the neutrino masses within our scenario, we need l/erl ~ IZ/xtI ~ 
\fe|J,\{.fn^i/'>^^■T)^ and a quite drastic finetuning between a and d in the mass matrix ©, 
in order to accommodate the solar mass-squared difference. Of course, the mass-squared 
difference AttIq ~ 10~^ eV^ of the quasi- vacuum oscillation solution [2] would require a 
stronger finetuning than AtUq ~ 10~^ eV^ of the large mixing angle solution. Examination 
of the Zee model with a general r2 and with small deviations of the neutrino mass matrix 
from the form (0) should reveal how much our finetunings could be relaxed. In the present 
work we have confined ourselves to prove that it is possible to reproduce the LMA MSW 
solution within the Zee model. 

In conclusion, even in the case that bimaximal mixing is ruled out, the Zee model will 
remain a viable and interesting scenario in order to accommodate the neutrino oscillation 
solutions of the solar and atmospheric neutrino problems. 
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A Charged scalar exchange and the 1-loop neutrino 
mass 

We proceed from the general Lagrangian 

- Cy{i^l, i, <^+) = {i^LARiR + 'M^Adi)^^ + h.c, (Al) 

where the Al^r are 3x3 coupling matrices in the case of 3 families and Il,r denote chiral 
charged fermion fields. The neutrino Majorana mass Lagrangian is defined by 

^.mass = ]^ylC-^M,UL + h.C. (A2) 

With the scalar interactions (jAlj) we obtain 

_y^i-ioop _ AlMjIIM^, Mf)A^ji + transp., (A3) 
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where we have defined 



1 



(47r)2 



e 



1 



+ ln(47r) + 7b - 1 + 



M2 -m' 



1 




) 



) 



(A4) 



We have used dimensional regularization; thus, e = (4 — n)/2 with n being the number of 
space-time dimensions, 'Je is Euler's constant and /i is an arbitrary mass scale. 
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